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ABSTRACT 

This paper examines tlie detailed sub-milliarcsecond polarization properties of an individual 
SiO maser feature displaying a rotation in polarization electric vector position angle of approx- 
imately ^ across the feature. Such rotations are a characteristic observational signature of cir- 
cumstellar SiO masers detected toward a number of late-type, evolved stars. We employ a new 
calibration method for accurate circular VLBI polarimetry at millimeter wavelengths, to present 
the detailed Stokes {/, Q, U, V} properties for this feature. We analyze the fractional linear and 
circular polarization as a function of projected angular distance across the extent of the feature, 
and compare these measurements against several theoretical models proposed for sharp rotations 
of electric vector position angle in polarized SiO maser emission. We find that the rotation is 
most likely caused by the angle between the line of sight and a projected magnetic field cross- 
ing the critical Van Vleck angle for maser propagation. The fractional linear polarization profile 
mi{6) is well-fit by standard models for polarized maser transport, but we find less agreement 
for the fractional circular polarization profile mc{0). 

Subject headings: masers - polarization - stars:individuaI(TX Cam) - starsmagnetic fields 



1. Introduction 

Circumstcllar SiO maser emission in commonly- 
observed transitions, such asi^e{l,2}, J — 1 — 0, 
orj/ = l, J = 2 — 1, is ubiquitous toward large- 
amplitude, lon g-period variable (LALPV) stars 
( Habin The excitation conditions for 

these masers place them in the near-circumstcllar 
environment (NCSE), between the stellar photo- 
spher e surface and the inner dust-formation ra- 
dius (|Elitzuij[l980h . but with the extent of the 
latter separation dependent o n stellar pulsation 
phase (jWittkowski et g/J 120071 ). The maser emis- 
sion is observed to be sign if icantly linearly po- 
larized (jTroland et all \l97^ iHerpin et ail l2006l ) 
but has a lower level of measured circular polar- 
ization (Barvain j g. Mcintosh, fc Predmord Il987 : 
Herpin et aLll200d ). The intrinsic high-brightness 
temperature and compact spatial st ructure of indi- 
vidual SiO maser components ( Moran et al. 19791 ) 



make them valuable scientific probes of the astro- 
physics of the NCSE. Very Long Baseline Interfer- 
ometric (VLBI) polarimetry, in particular, allows 
the properties of the NCSE to be imaged at sub- 
mas angular resolution by using the SiO maser 
components as tracers (jKemballl 120021) . Particu- 
larly important amongst the NCSE properties are 
the magnetic field magnitude and distribution, 
both at local and global stellar scales, and its as- 
sociated dynamical influence in this environment. 

Inference of magnetic field properties from the 
measured sub-mas SiO maser polarization however 
requires the solution of an inverse problem involv- 
ing the radiative transfer of SiO maser emission in 
full polarization. 

There remain differences in theoretical models 
for the transport of polarized maser emission in 
the limit of small Zeeman splitting, as applicable 
to the non-paramag netic SiO molecule (Elitzur 



2002; IWatsonI 12002). Both collisional 



1980 ) and radiative 



Elitzur 



Buiarrabal fc Nguven-Q-Rieu 



19811 ) mechanisms have also been proposed for SiO 
maser pumping. To complicate matters further, 
the physical conditions in the NCSE include com- 
plex phenomena such as shocks, anisotropic stellar 
illumination of maser cells, and likely both veloc- 
ity and magnetic field gradients, all of which affect 
the interpretation of maser polarization observa- 



tions. 

Spectral-line observations in full Stokes {/, Q, U, V} 
with sufficient angular resolution to resolve component- 
level SiO polarization properties can however pro- 
vide unique constraints on the theoretical and 
observational uncertainties enumerated above. In 
this paper, we report such measurements for an 
individual SiO maser feature displaying a distinc- 
tive observational signature, namely a rotation 
of approximately ^ in polarization electric vector 
position angle (EVPA) across the component at 
the inner shell boundary. 

The presence of such components was noted 
in the first VLB I polarimetry of SiO maser s to- 
ward TX Cam (jKemball fc Diamond! [l997t) but 
has also been noted in some subsequent obser- 
vations of other circumstellar SiO mase r sources 
(|Desmurs et aZ.ll2000l : ICotton et a/.ll2006h . 

Several distinct theoretical explanations have 

been p roposed for such features. lAsensio Ramos. Landi Degl'Innoceiiti 
(j2005l ) performed an analysis of the impact 
of radiative anisotropy on the excitation and 
linear polarization of unsaturated circumstel- 
lar SiO masers. In the non-magnetic case, 
linear polarization tangential to the projected 
shell is predi cted, as described also in ear- 
lier work by IWestern fc Watsoiil ( 19831 ). In 
evaluating possible causes of large EVPA ro- 
tations from this default tangential orientation. 
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(|2005l ) considered the impact of the Hanle effect in 
the presence of magnetic fields of order 10-100 mC 
They find that although the Hanle effect can pro- 
duce substantial EVPA rotations for a restricted 
range of magnetic field geometries, the required 
conditions produce a predicted suppression of the 
masing effect; it is therefore ruled out in their 
analysis. The authors conclude instead that the 
EVPA rotation is most likely caused by a change 
in the radiation anisotropy conditions across the 
maser component, from radial anisotropy to tan- 
gential anisotropy moving outward from the inner 
shell radius. In this model, the magnetic field is 
not dynamically significant. 

A second theoretical explanation involves the 
transport of linearly polarized maser emission in 
the presence of a magnetic field inclined at an an- 
gle to the line of sight. Foundational work on po- 
lar ized maser radiation transf e r was undertaken 
bv lColdreich. Keelev. fc Kwan ( 19731 ). In general. 



2 



the applicability of polarized maser transport so- 
lutions and their derivation is characterized by the 
relative magnitude of the maser stimulated emis- 
sion rate R, radiative decay rate F, Zeeman rate 
gfl, and bandwidth of the maser am plified radia- 
tion Aw (jGoldreich. Keelev. &: Kwa n 1973). For 
the case of saturated linear masers isotropically 
pumped over quantum number m, in the asymp- 
totic limit Aoj ^ gfl :s> R ^ T, these authors 
derived a solution in which the measured EVPA 
will rotate by ^ as the angle 9 between the mag- 
netic field direction and the line of sight crosses 
the critical Va n Vleck angle Op — 55° , defined 



as sin-^ (^f = j (|Goldreich. Keelev. &: Kwanlll973l) 



This change in 9 can therefore be gradual but still 
result in an EVPA rotation of ^. For 9 < Op the 
measured polarization EVPA is parallel to the pro- 
jected magneti c field direction, and perpendicu - 
lar for 9 >9f (|Goldreich. Keelev. fc KwaJ]ll973^ . 
lElitzuil (|2002[ ) noted that a transition across 9p 



could be responsible for the reported | rotations 
in EVPA across SiO maser features. We note 
that such an EVPA rotation has been detected in 
the polarized water maser emi ssion toward W43A 
( Vlemmings fc Diamondl 20061 ). 

A third theoretical explanation for significant 
EVPA rotation ac ross SiO ma ser components has 
been presented by lSoker ( 2002h . as part of a larger 
analysis of the likely dynamical influence of mag- 
netic fields in the evolutio n of late-type, evolved 
stars ( Soker fc Zoabi 20021) . These authors do not 
support the model of a globally-organized mag- 
netic field that is active during AGB evolution and 
acts dynamically to shape later planetary nebula 
geometry. However their work does support cool 
spots on AGB s tellar surfaces as sites o f enhanced 
dust formation ( Soker fc Clavton 19991) wit h pos- 
sible local magnetic fields of order 1-10 G ([Soker 
20021 ). The predicted local magnetic field geom- 
etry is radial above the c ool spot, bu t tangential 
closer to the photosphere (|Sokeijl2002l ). The NCSE 
is traversed by shocks ca used by the p ulsation of 
the central LALPV star (|Bowenlll988f ). The con- 
trast in magnetic field direction is expected to be 
enhanced by post-shock compression (.Soker 2002.) . 
The intrinsic change in magnetic field orientation 
from tangential to radial over a short projected 
angular distance near a cool photos phere spot is 
accordingly proposed bv lSokeii (|2002l) as the cause 
of EVPA changes of ^ across individual SiO maser 



components. 

In this paper, we utilize recent algorithmic de- 
velopments that allow accurate circular polariza- 
tion measurement in mi llimeter- wavelength spec- 
tral VLBI observations (jKemball fc Richteiil201l[ ) 
to analyze the linear and circular polarimetric pro- 
file across an individual SiO component exhibit- 
ing an EVPA rotation of approximately f . We 
measure the fractional linear and circular polar- 
ization profiles across the feature, and assess them 
against the theoretical models discussed above. 
We find that the EVPA rotation is most likely 
caused by the angle 9 between the line of sight 
and a magnetic field crossing the critical angle 9f 
noted above. 

The paper is organized as indicated below. The 
observations are described in Section 2 and their 
analysis and resulting science products presented 
in Section 3. The science results are discussed in 
Section 4, and conclusions presented in Section 5. 

2. Observations 

The current data are from a single epoch of 
a more extensive monitoring campaign using the 
Very Long Baseline Array (VLB A), operated by 
the NRACQ, to image thet/ = l, J:=l-0 cir- 
cumstellar SiO maser emission toward the Mira 
variable, TX Cam. These observations were 
conducted at sub-milliarcsecond angular resolu- 
tion and in full polarization Stokes {I,Q,U,V}, 
over several pulsation periods of the central star. 
The larger sur vey has been published in to- 
tal intensity by Diamond fc Keinballl ( 2003 ) and 



Gonidakis. Dianiond. fc Kemballl~2010l )^ and in 
linear polarization bv lKemball et al. Rioog). 

The current epoch was observed under VLBA 
project code BD46AQ. The observations were 
scheduled on 6 February 1999 (MJD 51215), from 
0'' UT to 8'' UT using all ten VLBA antennas 
plus an additional single VLA antenna. A total 
of 6.5 hours were assigned to the 43 GHz obser- 
vations, divided into on-source scans of 13-minute 
duration, and comprising the following total in- 
tegration times: i) seventeen scans on the target 
source, TX Cam; ii) seven scans on the contin- 
uum extragalactic calibrator J0359-I-509; and iii) 



^The National Radio Astronomy Observatory is a facility of 
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one scan each on the continuum extragalactic cah- 
brators 3C454.3 and J0609-157. The scans on TX 
Cam and J0359+509 were distributed as evenly 
as possible to maximize uw— coverage. Approx- 
imately 86.7% of the elapsed schedule time was 
on-source, the remainder was allocated for an- 
tenna slewing, system initialization, or to balance 
the tape resources allocated at that time. 

The V = 1,J = 1 — SiO maser transition 
was observed in a 4 MHz baseband, and centered 
on a systemic LSR velocity of -1-9 kms~^ for TX 
Cam. The adopted rest frequency was 43.122027 
GHz. No real-time Doppler tracking was employed 
about the mean Doppler shift computed for the 
mid-point of the schedule and the array; these 
corrections were applied in post-processing, con- 
sistent with sta ndard spectral-line VLBI practice 
(Diamond Il989l) . The data were sampled in 1-bit 
quantization, and correlated in full polarization 
over a maximum possible 128 frequency channels 
available at that time, yielding a nominal channel 
width of 31.25 kHz. The correlator accumulation 
interval was 4.98 s. 

3. Results 

The da ta were reduced using the method de- 
scribed bv lKemball fc Richter (2011) for accurate 
circular VLBI polarimetry at millimeter wave- 
lengths. The primary science products resulting 
from the reduction are image cubes in each of 
Stokes {I,Q,U,V}, at a pixel spacing of 50 ^as 
(2048 X 2048) on the image tangent plane, with 
one image per sampled frequency channel over the 
inner 113 frequency channels in the spectrum. A 
common restoring beam of size 540 x 420 //as at a 
position angle of 20° was adopted a cross all epochs 
of the larger survey, as described bv lKemball et al. 



()2009l) . The absolute electric- vector position angle 
(EVPA) of the linearly-polarized emission was es- 
tablished from associated VLA observations, rel- 
ative to the primary polar ization EVPA c alibra - 
tor, 3C138, as described bv lKemball et al\ (|2009f ). 



The residual error in the absolute EVPA determi- 
nation is estima ted to be ~ 10° — 20° peak-to-peak 
20091 ). 



(Kemball et al 



The zeroth moment, over frequency, of the full 
Stokes I image cube is shown in Figure [1] The 
counterpart Stokes V image is shown in Figure 2. 
The linearly-polarized intensity P — \/ [Q^ + t/^). 



derived directly from the zeroth-moment images 
in Stokes Q and [/, is depicted in Figure 3. In 
Figure 4, the Stokes / image is shown overlaid with 
vectors proportional in length to P and drawn at a 
position angle equal to the absolute EVPA of the 
linearly-polarized emission. 

This paper concerns a feature located near 
relative coordinates (-15, -7.5) mas in Figure 4 
that has an approximate ^ rotation in EVPA 
across the feature. We note that the image cubes 
lack absolute astrometric coordinates due to the 
use of VLBI phase self-calibration. This compo- 
nent, near the southwest circumstellar boundary, 
is shown in Stokes / and P in Figure El averaged 
over frequency as in earlier Figures. As noted ear- 
lier, it is not uncommon to find adjacent individual 
SiO components with large (or perpendicular) ro- 
tations in relative EVPA. In the current epoch, we 
choose this particular individual feature for fur- 
ther analysis because it spans a relatively large 
range in velocity and has good SNR. 

The individual frequency channel images across 
this feature, labeled by line-of-sight velocity (in 
km/s), are plotted in Stokes / and P in Figure [HI 
and in Stokes V in Figure [71 

In the absence of knowledge of the astromet- 
ric position of the central star relative to the 
SiO emission, we adopt the approximation that 
radially-extended maser features point back to the 
photosph ere, an assumption for SiO masers dis- 
cussed by [Zh^g^^el^ (2011). The corresponding 
vector obtained by fitting a straight line to the 
projected coordinates of the peak Stokes / com- 
ponent brightness across the velocity extent of the 
^ EVPA rotation feature is drawn as an arrow at 
upper left in Figures jSJ [51 and [71 and annotated 
accordingly as the approximate assumed direction 
toward the photosphere. 

The integrated mean-intensity spectrum, com- 
puted across the full image region in Figure |S] 
enclosing the | rotation feature, is plotted in 
Stokes /, Stokes V , and linearly-polarized inten- 
sity P = ^JQ^^^^ in Figure HI 

4. Discussion 

The large-scale morphological properties of the 
v — 1, J — 1 — SiO maser emission distribution 
at this epoch in Stokes / and P are consistent with 
the summary results for the broader monitoring 
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campa ign, as discussed b y iDiamond fc Kemball 
(|2003l ) and iKemball et aZI ([20091). The total in- 
tensity distribution, depicted as an average over 
frequency in Figure [TJ shows the projected ring- 
hke shell morphology frequently found for cir 



cumstellar SiO maser emission (jPiamond et al. 
1994I ). The corresponding linearly-polarized in- 
tensity, shown in Figures 3 and 4, has the charac- 
teristic tangential distribution of EVPA found to 
be per sistent across the br oader monitoring cam- 
paign ([Kemball et al. [ l2009l) . However this tangen- 
tial polarization morphology is not universal for 



circu mstellar SiO masers as a class ([Cotton et al. 
20091 ). 



The individual component studied in this paper 
is located at the projected shell boundary in the 
southwest region of the overall SiO maser distri- 
bution. As shown in Figure [5l there is a rotation 
of approximately ^ in EVPA across the compo- 
nent, with a tangential orientation on the inner 
shell boundary and a radial orientation at a larger 
projected angular distance from the central star 
(which lies toward the northeast in Figure [5]) . 

The component is elongated in total intensity 
along the radial axis. This is consistent with 
tangential amplification o f the underlying maser 
emission, as described bv [Diamond et al. ( 1994 ). 
and consistent with radial shock acceleration pre- 
dicted in the circumstellar shells of LALP V stars 
([Humphrevs et a/.|[2002[: [Crav et a/.|[2009[ ). The 
total intensity contour images of each frequency 
channel across the feature (Figure[6|) show that the 
projected center of the Stokes / emission moves in- 
ward toward the central star with decreasing LSR 
velocity, as expected for a positive radial gradient 
in the outfiow velocity at the position of this SiO 
maser component in this shell. 

The linear polarization channel images shown 
in Figure |6] span the feature in velocity. They in- 
dicate an abrupt transition in EVPA of ^ near 
Vlsb. ^ +6.2 km/s and a corresponding local 
minimum in linearly-polarized intensity near this 
point. The corresponding circular polarization 
channel images in Figure [7] show a sharp decline in 
peak Stokes V with decreasing LSR velocity, i.e. 
in the direction toward the inner projected shell 
boundary. 

The fractional linear (m;) and circular (rric) po- 
larization magnitude profiles across the feature are 
plotted in Figure [9l for all measurements with an 



SNR exceeding 3. The magnitudes of the frac- 
tional polarizations are measured at the single 
pixel position of maximum Stokes / in each chan- 
nel image. The a;— ordinate in this plot is the pro- 
jected angular separation d (mas) from the compo- 
nent peak in the channel image at Vlsb. — +7.91 
km/s in the upper left panel of Figure El In this 
sense, the projected angular separation increases 
toward the central star - with decreasing Vlsr. due 
to the radial velocity gradient discussed above. 

The single-pixel measurements of to/ include 
the effects of spatial linear depolarization arising 
from convolution by the synthesized beam during 
image formation. The magnitude of this effect can 
be assessed for a given deconvolved source compo- 
nent size am, a synthesized geometric beamwidth 
cTft, and an adopted linear rate of change of EVPA 
a = X with angular spatial scale in the image. 
The angular sizes a are expressed here as the 
full width at half maximum intensity (FWHM). 
The components in Figure [5l with the largest ap- 
parent values of a are at Vlsr — +6.18 km/s 
and Vlsa — +5.96 km/s. The mean decon- 
volved source component size across the minor axis 
(which is almost perpendicular to the EVPA), for 
these velocities is am ~ 1-2 mas. The geomet- 
ric synthesized beamwidth for the current data is 
cTfc ~ 0.48 mas (see above). The linear beam de- 
polarization arising in this case is approximatly 
toJ ~ ^^IJ^f^™' ~ with a — Adopt- 

ing — f and Ax^ — j produces /? = 0.94 

and (3 = 0.98 respectively, with values closer to 
unity for the other components in Figure [HI due 
to significantly lower apparent values of a. Thus, 
for the single-pixel measurements of to; presented 
here, spatial beam depolarization is not believed 
to substantially affect the current analysis. 

For the case of the single-pixel fractional circu- 
lar polarization measurements, TOc, it is similarly 
appropriate to consider the effect of depolarization 
arising from averaging over frequency. As noted 
earlier, the nominal channel increment in the data 
is 31.25 kHz; in uniform spectral weighting the ef- 
fective FWHM spectral resolution is a.^ ~ 0.26 
km/s. The mean observed FWHM component 
line-width in Stokes / over the feature is aots ^ 0.9 
km/s. To first-order, for a classical Zeeman "S- 
curve" the peak Stokes V obeys the proportion- 
ality Vmax oc for a Zeeman splitting 
and component line-width a. Convolution by the 
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instrumental frequency response therefore scales 
V„iax by a factor -q ^ y^l — (^^)^- For the val- 
ues of (Tt, and CTo&s for the current data, 77 ~ 0.96, 
which does not significantly affect the results pre- 
sented in the current work. This calculation is 
illustrative only however, as the Stokes V com- 
ponent profiles do not take simple Zeeman form, 
as evident from the integrated spectrum plotted 
in Figure HI However, we believe the calculation 
is nonetheless representative of the magnitude of 
the circular depolarization effect. We discuss the 
broader question of depolarization caused by line- 
of-sight integration below. 

The asymptotic linear polarization solution for 
saturate d, m— isotropic line a r mas ers derived by 
[Goldrcic h. Keelev. fc Kwaiil (1973) for the case 
where Aw 3> gfl ^ R ^ T predicts a fractional 
linear polarization dependence on the angle 9 be- 
tween the ma gnetic field and the line of sight ta k- 
ing the form (iGoldreich. Keelev. fc K-JMlil973l) : 



miie) 

mi{9) 



3 sin^ 6 

for e> Ob 

If — 



3 sin^ e 

1 for e <eB 



(1) 



where tan^^s = ^. In this solution the measured 
EVPA rotates by § at 6>_f, wher e sin^6'F = | 
(|Goldreich. Keelev. \ KwanI Il973l ). We denote 
this solution as GKK in what follows. We note 
that the GKK solution was derived specifically for 
a J = 1 — transition, as applies to the transition 
observed here. 

We model the fractional linear polarization data 
across the SiO maser feature by fitting the lowest- 
order polynomial form of 0(d) on the projected 
angular distance d across the feature that yields 
a reasonable fit. For the current feature this is a 
quadratic fit: 



(2) 



where dj — 2.822 mas is the projected angular 
separation at which the EVPA rotates by ^ in the 
measured data - here chosen to be the channel at 
yhSR = 6.18 km/s. The free parameters in the 
fit are a and b. In Figure [9] we plot the best chi- 
square fit of the measured nii data to the GKK 
solution as a dashed line. The associated solution 
for 9{d) is plotted in Figure [TOl 



The fit to mi shows broad agreement with the 
functional form of the GKK solution. In this 
model there is a gradual change in the angle 9 be- 
tween the line of sight and the projected magnetic 
field over the range 37.5° to 70° across the feature 
(with increasing d). The angle 9 crosses the Van 
Vleck angle 9f a.t d — df near Vlsr — +6.2 km/s. 
We note that for this feature the lower bound of 
37.5° is consistent with t he lack of a sta ble solution 
for 9 < 9b predicted bv lElitzuil (| 19961) . 

If we assume that m; {9) originates from a GKK 
solution then we can derive the dependence of frac- 
tional circular polarization on 9 by direct inversion 
of Equation [1] as: 



cos 9 = . / 1 — 



3(mi + l) 



(3) 



We plot the resulting dependence of measured 
fractional circular polarization on 9 in the form 
mc(cos6') in Figure[TTJ Here, the measured values 
of mi are used to compute cos 6* using Equation [3j 
the matching values of mc are then plotted against 
the derived cos 9 abscissa. The plot shows an in- 
crease of TOc with cos 9 over a relatively narrow 
range of angles 9 for which the projected magnetic 
field is closest to the line of sight. 

We can assess the measured component-level 
polarization properties presented here relative to 
the different theoretical models that have been 
proposed to explain EVPA rotations of approx- 
imately 5 across individual SiO maser features, 
described in the Introduction. However, there are 
important caveats that apply; we know this region 
has magnetic field and velocity gradients and is in 
a shock-traversal region; many theories were de- 
veloped within more idealized conditions, by ne- 
cessity. Furthermore, our observations are inte- 
grated over the synthesized angular beamwidth 
along lines of sight within this complex environ- 
ment . We do not believe this affects substantially 
the primary conclusions of this paper however, but 
more complex three-dimensional modeling of this 
region is planned in future work. 

Th e non-magnetic {gf!, — 0) model of lAsensio Ramos. Landi Degl'In 
( 20051 ) explains the EVPA rotation as due to 
changes in radiative isotropy across the SiO maser 
feature. We believe that the current data do not 
provide support for this model. If the tangen- 
tial linear polarization at the inner edge of the 
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shell is predominantly caused by m-anisotropic 
pumping in a non- magnetic environment, then 
the measured circular polarizati on would hav e to 
arise from non-Zeeman effects ( Watson 2009() : in 
this case the inter-conversion of linear polariza- 
tion to circular polarization is due to a change 
in optical axes along the propagation path. In a 
non-magnetic model, we might expect from first 
principles that non-Zeeman circular polarization 
would be found preferentially at positions with the 
greatest EVPA rotation rate; however the circular 
polarization is not pronounced at this position in 
our data. We also note that the greatest frac- 
tional linear polarization is measured in our data 
at the feature position that is furthest from the 
central star. Further, the relatively close agree- 
ment with the GKK functional form for mi{0) 
argues against a model of EVPA rotation due pri- 
marily to changes in radiation isotropy. For these 
reasons we believe that our data do not provide 
support for the EVPA rotation model described by 
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(120051). 

In considering models with non-zero magnetic 
fields g^l 7^ 0, we believe that contemporary ob- 
servational evidence suggests that circumstellar 
SiO masers are in the regime qil ^ R ^ F 



gfl > R :$> T (|Kemball et all l2009l ). Evi 



dence in support of the partial subsidiary con- 
dition gfl R is provided by Watson (2009). 
This condition excludes intensity-dependent non- 
Zeeman circular polarization that is possible if 
gVL ^ R (Ncdoluha & Watson 1994). The condi- 
tion gVl 3> R also ensures that the magnetic field 
is always either parallel or perpendicular to the 
measured EVPA (as in the GKK model presented 
above), but not necessaril y at the same value of 
the Van Vleck angle Of dWatsonl [2009h or with 
the same functional form for mi{9). 

Other sources of non-Zccman circular polariza- 
tion are possible for maser models with gfi ^ 0, 
due also to differences between the direction of 
maser linear polarization and the optical axes 
along the propagation path, in this case caused 
by chang es in magnetic field direction or Farada y 
rotation (jWiebe fc WatsonI Il998l : IWatsonI l2009l ). 
It is important to consider whether these sources 
of non-Zeeman circular polarization could explain 
our current observations. For this mechanism, 
we expect a correlation between fractional circu- 



lar polarization and fractional linear polarization, 
and that the greatest circular polarization will oc- 
cur at positions with the gr eatest rate of change 
of magnetic field direction (Watson 2009[ ). Our 
current data are inconclusive on this point. There 

is a partial correlation between rric and mi, and 

2 

within the constraint rr ic < put forward by 
Wiebe fc WatsonI (|l998l ) for this mechanism, but 



over too small a number of samples to allow a 
robust statistical conclusion from these data (see 
Figure [9] for reference). In contrast, the magnetic 
field angle gradient plotted in Figure [TU] (derived 
from the GKK fit) is anti-correlated with mc. 

The linear polarization profile shown in Fig- 
ure [9] can also be examined in terms of the infor- 
mation it provides on predicted maser saturation 
in different models of maser polarization propaga- 
tion. The GKK solutions assume st rong satura- 



tion R^T. In the work of lElitzuij (|1996f ) a linear 
polarization solution of GKK form can be attained 
well before saturation; in contrast the models of 
Watson k. Wvldl (|200l[ ) require strong saturation 
in a uni-directional J = 1 — maser to achieve the 
value of TO; 0.7 reported here. However, we do 
note that the latter model assumes m— isotropic 
pumping. In light of these aggregate predictions, 
we believe our linear polarization profile is consis- 
tent with these masers being saturated. 

This admits an interpretation in which our lin- 
ear polarization data are explained as resulting 
from saturated maser emission i? 3> F from a re- 
gion threaded by a magnetic field that changes di- 
rection smoothly relative to the line of sight across 
the maser feature. The angle between the mag- 
netic field and line of sight is approximately 37.5° 
at the furthest point from the star, and approx- 
imately 70° at the inner shell boundary of the 
feature. The relative orientation crosses the criti- 
cal Van Vleck angle Op within this region, causing 
an ab rupt § change in measured EVPA dElitzur 
2002) . We know from earlier work (Kemball et al. 



l200ilCotton. Perrin. fc LopedlioOSt ) that Individ 
ual maser motions appear to be influenced by in- 
dividual magnetic field lines. There are also very 
plausible mechanisms for producing such field cur- 
vature, even if only local fields are considered such 
as those associat ed with prop osed AGB cool spots, 
as noted earlier ( Sokeiil2002l) . 

Our measured circular polarization profile 
raises several theoretical issues. The GKK solu- 
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tion used earlier makes no prediction about circu- 
lar polarization - it is identically zero at line center 
in their asymptotic solutions (Watson 2009). The 
foundational work of GKK has been generalized 
to a wider range o f pararneter space in subsequent 



studies however (lElitzuij |1996t IWatson fc Wyld 



200 II ). and we examine those predictions here. 

The current data show a possible increase in 
rric with cos 6, perhaps of linear form, but with 
a sharper than expected fall-off with increasing 
9. These conclusions are tempered by the lim- 
ited number of points however, and clearly fur- 
ther data are needed. A cos 61 functional depen- 
dence occurs for thermal Zeeman emission, and for 
the case of widely-separated Zeeman mase r com- 
ponents (iGoldreich. Keelev. fc Kwan 1973 ). The 
work of lElitzuij (|l996l ) predicts a rric oc cos 9 
dependence for t he sm all Zeeman-splitting case. 
Watson fc Wvldl (|200l[ ) predict a circular polar- 
ization profile over 9 that changes in shape as a 
function of degree of saturation ^ ■ A linear rela- 
tion between rric and cos 6 is not predicted in this 
theory except for highly unsaturated emission; for 
increasing saturation uric increases sharply toward 
a peak in the region cos 6 < .5, then declines 
to ze ro as expected a,t 6 = j (jWatson fc Wyld 



20011 ). However, as noted earlier, this model is for 
m— isotropic pumping. Further observations are 
needed to elucidate the nature of the circular po- 
larization profile. 

5. Conclusions 

We have examined an individual SiO maser fea- 
ture that shows an EVPA rotation of approxi- 
mately ^ over the projected extent of the emis- 
sion. From our analysis, we find: 



1. 



The fractional linear polarization, as a func- 
tion of the angle 9 between the line of 
sight and the magnetic field, reproduces 
the functional form of the maser polariza- 
tion radiative transfer solu t ion d erived by 
Goldreich. Keelev. fc K^anl (|l973f ) for the 



case Aw ^ g^l R ^ T. A simple model 
is adopted for the relation between 9 and 
the projected angular separation across the 
feature. 

2. The rotation in EVPA by approximately ^ 
is explained within this model by 9 crossing 



the Van Vleck angle 9f near the mid-point 
of the component. 

3. Further observations are needed to clarify 
the functional dependence of circular polar- 
ization on cos 9. 



We are grateful to our colleagues for their com- 
ments on earlier drafts of this paper. We also 
thank the journal referee for comments that im- 
proved the clarity and scientific content of the 
manuscript. 

Facilities: VLBA. 
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Fig. 1. — Stokes / contour image of the v = 1, J = 
1 — SiO maser emission toward TX Cam, plotted 
as the zeroth moment over frequency of the image 
cube. In these averaged units, the contour levels 
are at {-10, -5, 5, 10, 20, 40, 80, 160, 320} x a, 
where a is the off-source rms of 2.1982 mJy/beam. 
The angular coordinates are in mas from the cen- 
ter of the sub-image enclosing the projected SiO 
maser ring. 
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Fig. 2. — Stokes V contour image of the 
u — 1, J = l — SiO maser emission to- 
ward TX Cam, plotted as the zeroth mo- 
ment over frequency of the image cube. In 
these averaged units, the contour levels are at 
{-160, -80, -40, -20, -10, -5, 5, 10, 20, 40, 80, 160}x 
a, where a is the off-source rms of 1.7113 
mJy/beam. The angular coordinates are in mas 
from the center of the sub-image enclosing the 
projected SiO maser ring, aligned with Figure [T] 
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Fig. 6. — Stokes / channel images, labeled with LSR velocity (km/s), plotted at contour levels 
{ — 12, —6, —3, 3, 6, 12, 24, 48, 96, 192, 384, 768}a, where a — 15.7 mjy/beam is the off-source rms in an early frequency 
channel in this sequence. Overlaid linear polarization vectors are drawn as in Figure [5] but with scale such that P= 
1.05 Jy/beam has length 1 mas. 
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Fig. 7. — Stokes V channel images, matcliing Figure (6] plotted at contour levels 

{-96, -48, -24, -12, -6, -3, 3, 6, 12, 24, 48, 96} x a, where a = 15.7 mjy/beam. 
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Fig. 3. — Stokes P contour image of the v = 
1, J = 1 — SiO maser emission toward TX Cam, 
plotted as the zeroth moment over frequency of the 
image cube. In these averaged units, the contour 
levels are at {7.5, 15, 30, 60, 120, 240, 480} x cr, 
where a is the off-source rms of 1.1954 mJy/beam. 
The angular coordinates are in mas from the cen- 
ter of the sub-image enclosing the projected SiO 
maser ring, aligned with Figure [T] 
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Fig. 4. — Stokes / zeroth-moment image over fre- 
quency, plotted at the contour levels in Figure [T] 
the overlaid vectors are drawn with position an- 
gle equal to the absolute EVPA of the underlying 
zeroth-moment linearly-polarized intensity P and 
with length proportional to P such that P= 16 
mJy/beam has length 1 mas. 
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Fig. 5. — A sub-image from Figure 4 enclosing 
the south-west component on the projected shell 
boundary with a ^ rotation in EVPA. The Stokes 
/ contour levels are as used in Figure 4, but the 
linearly-polarized zeroth-moment intensity P is 
drawn with vector length on a scale where P = 
26.7 mJy/beam has length 1 mas. 
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Fig. 8. — Mean-intensity spectrum in Stokes /, 
V, and P = y^Q"^ + over the full image region 
in Figure [5l enclosing the feature showing the | 
rotation in EVPA. This spectrum was obtained 
by integration over the associated interferometric 
image cubes. The Stokes V profile is magnified by 
a factor of three for clarity of presentation. 
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Fig. 10. — Solution for 6{d) {left axis) derived in 
the fit to the GKK model shown in Figure |9l The 
gradient of 9{d) with respect to d is plotted as a 
dashed line (right axis). 
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Fig. 9. — Magnitude of fractional linear and circu- 
lar polarization at the position of the peak Stokes 
I intensity in each frequency channel spanning the 
component shown in Figure [5l for all components 
with an SNR above three. The projected angu- 
lar distance d is measured as the linear separa- 
tion from the peak of the upper left channel im- 
a ge in Figure [HI The linear polarization model 



a ge m f igure lb) i ne nnear polari zation moaei 
of boldreich. Keelev. fc k"^(Ii973I) . fitted as de- 
scribed in the text, is plotted as a dashed line. The 
associated 9{d) relation is plotted on the upper 
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Fig. 11. — Fractional circular polarization as a 
function of cos 6, where 9 is the angle between 
the magnetic f ield and line of sight, compu t ed ac - 
cording to the Goldreich. Keelev. fc Kwan (197^ 

model, as cos 6* , /l — 3(„^^^;^) ■ 



a;— axis. 
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